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Semiconductor Device Simulation Method and Simxilator 

The present invention relates to an automated simulation method for 
determining the enhanced generation-recombination rate due to trap-to-band 
tunnelling in a semiconductor device using the Dirac coulombic tunnelling integral 
and to a simulator for carrying out the method. 

The enhanced generation-recombination rate due to trap-to-band tunnelling in 
a semiconductor device is considered to be responsible for a number of important 
effects in semiconductor devices, including the anomalous leakage current in 
polysilicon Thin Film Transistors (TFTs). Polysilicon TFTs are widely used for 
example in active matrix display devices. The anomalous leakage ciirrent of the 
TFTs can severely degrade the pixel voltage in such display devices. Thus, this is 
one example of the commercial importance of an automated simulation method and 
simulator of the type provided by the present invention. 

Having regard to the complexity and cost of the fabrication processes for 
manufacmring semiconductor devices, it is highly desirable if not essential for the 
design and performance evaluation of such devices to be undertaken using 
mathematical simulations, often referred to as modelling. It is clearly crucial for 
such modelling to be able to provide accurate calculation of the enhanced 
generation-recombination rate due to trap-to-band tunnelling, and hence the leakage 
current, within a semiconductor device. Consequently, considerable effort has 
previously been spent in developing methods of calculating the enhanced 
generation-recombination rate due to trap-to-band tunnelling in a semiconductor 
device. Such methods are embodied in computer programs which are sold as staple 
commercial products by or on behalf of their developers to designers and 
manufacturers of semiconductor devices. 

In a paper submitted in 1996 and published in 1997 (Solid State Electronics 
Vol.41, No.4, pp 575-583 1997) the inventors hereof presented a generation- 
recombination model for device simulation including the Poole-Frenkel effect and 
phonon-assisted tunnelling. The model is conveniently referred to as the Dirac 
Coulombic Tunnelling Integral, which is recited as equation 1 in figure 9 hereof. 
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As noted in the above mentioned published paper, the Dirac Coulombic 
Tunnelling Integral is applicable to semiconductor devices generally. However, as 
also noted above, an important category of semiconductor devices is thin film 
transistors (TFTs) and such a device will herein after be used for ease of reference, 
but by way of one example only of a semiconductor device. Similarly, for ease of 
reference, the leakage current in a TFT will be referred to herein as a non-limiting 
example of the enhanced generation-recombination rate due to trap-to-band 
tunnelling in a semiconductor device. 

Figure 1 hereof is a graph showing voltage current characteristics of a 
polysilicon TFT. As seen in Fig. 1 when V^s is high (5.1V), the leakage current 
(lus) increases with decreasing V^s (below OV). The magnitude of this leakage 
current poses as a significant problem, for example, when the TFT is employed as a 
switching pixel transistor in active matrix LCDs. Several field-assisted generation 
mechanisms have been proposed to explain this 'off' current, 

A quantitative analysis of the leakage current in polysilicon TFTs based on 
the combination of current-voltage measurements £ls a function of temperature and 
2-D simulations was carried out as long ago as 1995, This analysis shows that the 
dominant generation mechanism is pure trap-to-band tunnelling below 240K and 
phonon-assisted trap-to-band tunnelling at higher temperamres. The need to include 
Poole-Frenkel (PF) barrier lowering in trap-to-band phonon-assisted tunnelling was 
demonstrated for polysilicon pn junctions as long ago as 1982. Apart from further 
enhancing the emission rate for trap-to-band phonon-assisted tunnelling, the PF 
effect also plays a significant role in enhancing pure thermal emissions at low 
fields. 

The PF effect consists of the lowering of a coulombic potential barrier due 
to the electric field applied to a semiconductor. For a trap to experience the effect, 
it must be neutral when filled (charged when empty). Such a trap potential is long 
ranged and is often referred to as a coulombic well. A trap that is neutral when 
empty will not experience the effect because of the absence of the coulomb 
potential. Such a trap potential is short ranged and is known as a Dirac well. 
Without the PF effect, the calculated emission rate is at least one order of 
magnitude lower than what is needed to fit the experimental data in polysilicon 
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TFTs. By iising a 2-D simulator, based on the known trap-to-band phonon-assisted 
tunnelling model (Hurkx et a/) available since 1992, the inventors hereof were 
unable to simulate the leakage currents in polysilicon TFTs accurately. This is 
because the conventional model takes only into accoimt Dirac wells and deliberately 
neglects the PF effect. Moreover, the many attempts at modelling of trap-to-band 
phonon-assisted tunnelling inclusive of the PF effect, which have been made since 
development of the original theory by Vincent et al in 1979, do not address a key 
problem: the implementation in a device simulator. The work by Vincent et al in 
1979 gives evidence both theoretically and experimentally that the electric field in a 
junction has a large influence on the thermal emission rate of deep levels (mid-gap 
states). This influence can be quantitatively explained in a model of phonon- 
assisted tunnelling emission. Tunnelling is very sensitive to the barrier height and 
is therefore expected to be considerably affected by the PF barrier lowering. 

As noted above, in 1996 the present inventors presented a new quantum 
mechanical tunnelling generation-recombination (G-R) model, conveniently referred 
to as the Dirac Coulombic Tunnelling Integral, which takes into full rigorous 
account the PF barrier lowering and is suitable for implementation in a device 
simulator. This G-R model is consistently formulated for the entire range of 
electric fields and temperatures. At high fields, the dominant mechanism is found 
to be trap-to-band phonon-assisted tunnelling inclusive of the PF effect; while at 
low fields, the model will reduce to that of the standard Shockley-Read-Hall (SRH) 
thermal G-R. 

However, a practical implementation of the model into commercial device 
simulators has not imtil this invention been presented despite the length of time 
since presentation of the Dirac Coulombic Tunnelling Integral and despite the high 
commercial value of such a practical implementation. 

According to a first aspect of the present invention there is provided an 
automated simulation method for determining enhanced generation recombination 
rate due to trap-to-band tunnelling in a semiconductor device using the Dirac 
coulombic tunnelling integral, comprising the steps of: 

assigning the variable C to the ratio of the Poole-Frenkel barrier lowering 
energy (AE^) divided by the energy range for which tunnelling can occur (AEJ; 
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assigning the value (C+ 1)/2 to a variable v and performing a second order 
Taylor's series expansion of the Dirac coulombic tunnelling integral around v to 
determine a maximum value (u^aJ for the variable u of the integral; 

determining if the value for u^j^x is less than C, is between C and 1 or is 
more than 1 ; 

assigning the value of C to the variable v if is less than C; 

assigning the value of u^ to the variable v if is between C and 1 ; 

assigning the value of 1 to the variable v if u^ is more than 1; 

reducing the Taylor's series expansion of the Dirac coulombic tunnelling 
integral to an error function; 

reducing the error function to simple exponential functions by applying 
rational approximations to the error function; and 

calculating the enhanced generation recombination rate due to trap-to-band 
tunnelling in a semiconductor device using the said simple exponential functions. 

According to a second aspect of the present invention there is provided a 
simulator for determining enhanced generation recombination rate due to trap-to- 
band mnnelling in a semiconductor device using the Dirac coulombic tunnelling 
integral, comprising: 

means storing a variable C having a value equal to the ratio of the Poole- 
Frenkel barrier lowering energy (AE^) divided by the energy range for which 
tunnelling can occur (AEJ; 

means which assign the value (C + 1)/2 to a variable v and perform a second 
order Taylor's series expansion of the Dirac coulombic tunnelling integral around v 
to determine a maximum value (u^) for the variable u of the integral; 

means which determine if the value for u^ is less than C, is between C and 
1 or is more than 1 ; 

means which assign the value of C to the variable v if is less than C; 
means which assign the value of u^ to variable v if u^ is between C 

and 1; 

means which assign the value of 1 to the variable v if u^^ is more than 1; 
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means storing simple exponential functions derived from applying rational 
approximations to an error function obtained by reducing the Taylor's series 
expansion of the Dirac coulombic tunnelling integral; and 

means which calculate the enhanced generation recombination rate 
due to trap-to-band tunnelling in a semiconductor device using the said simple 
exponential functions. 

Although possible in theory, even with the fastest computers available today, 
using niunerical integration methods to solve the Dirac coulombic tunnelling 
integral for every trap level in every element in a finite element package will take 
so much computation time that such an approach can not be used for a commercial 
and practicable implementation. The present invention enables such a conmiercial 
and practicable implementation. Moreover, the present invention can provide high 
levels of accuracy in the automated modelling of the enhanced generation- 
recombination rate due to trap-to-band tunnelling in a semiconductor device. 

Embodiments of the present invention will now be described in more detail 
by way of example only and with reference to the accompanying drawings, in 
which: 

Figure 1 is a graph illustrating leakage current in a TFT; 

Figure 2 shows two graphs illustrating a low field condition; 

Figure 3 shows two graphs illustrating a moderate field condition; 

Figure 4 shows two graphs illustrating a high field condition; 

Figure 5 is a graph illustrating the performance of the present invention; 

Figure 6 is a graph illustrating the performance of the present invention; 

Figure 7 is a graph illustrating the performance of the present invention; 

Figure 8 is a graph illustrating the performance of the present invention; and 

Figure 9 lists equations useful in explaining embodiments of the present 

invention. 

Starting from the Dirac coulombic tunnelling integral, given as Eqn.(l) in 
figure 9, it is possible for the integral to be expressed in the form shown as 
Eqn.(2). Applying the simplification shown in Eqn. (3) leads to Eqn. (4). Eqn.(4) 
can be further simplified to give Eqn. (5), where the function /(u) is as defined in 
Eqn.(6). 
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Most of the contribution to the expression of Eqn.(5) can be expected to 
arise from the condition when the function ^u) is largest. This condition should be 
considered for the three cases which have distinct characteristics namely. Case 1 
when the electric field is low; Case 2 when the electric field is moderate and Case 
3 when then the electric field is high. These three distinct conditions are illustrated 
in figures 2, 3 and 4 respectively. 

Figure 2 is based on a low field value (F) of F = IxlO^Vm"' foj AE„ = 
0.5eV. The maximxrai of/(u) occurs at C, giving the maximum contribution to 
^IfMJ at C as well. 

Figure 3 is based on a moderate field value (F) of F = VxlO^Vm'* for AE^ 
= 0.5eV. The maximum of/(u) occurs between C and 1, giving the maximum 
contribution to exp[f(u)] between C and 1 as well. 

Figure 4 is based on a high field value (F) of F = IxlO^Vm"^ for AE^ = 
0.5eV. The maximum of/(u) occurs at 1, giving the maximum contribution to 
^IfM] at 1 as well. 

In a practical implementation it is convenient to consider Case 2 first. That 
is to determine the u^ for maximtun/(u) between C and 1. If u^ is less than C, 
then the condition of Case 1 is met. If u^ is more than 1 then the condition of 
Case 3 is met. Using a Taylor's series expansion, ^u) can be approximated by a 
second-order series expansion around v, where v = (C4- 1)/2 as a reasonable 
estimate of where u^ is likely to occur. Thus is derived Eqn,(7), assimiing^v) is 
as set out in Eqn.(8) and/(v) and f\v) are as set out in Eqn.(9) and Eqn.(lO) 
respectively. Further, this enables y(u) to be rewritten as in Eqn.(ll). 

From Eqn.(ll),^u) can be differentiated once and the expression equated 
to zero, to obtain a stationary point, thus following the Eqn.s(12) to obtain a value 
for Ujj^. If this value of u^ is less than C then Case 1 exists. Then it is 
permissible to set v = C since that is where the maximum ^(u) occurs. Then a 
second order Taylor's series expansion is performed about v = C. All of the 
necessary equations have been established, so it is only necessary to set v = C in 
the simulator when u^ is less than C, so as to obtain an approximate second order 
expansion offiyi) about C. 
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If Ujna, is between C and 1, Case 2 exists. In this condition, it is only 
necessary to set^v = in the simulator, so as to obtain an approximate second 
order expansion of^u) about Hj^. Here, of course, v = u^ is no longer assumed 
to be equal to (C + 1)/2. 

If u^ is more than 1, Case 3 exists. In this condition, it is only necessary 
to set V = 1 in the simulator, so as to obtain an approximate second order 
expansion o{f{u) about 1, 

With the appropriate value of v determined, Eqn.(ll) is simplified to lead to 
Eqn.(13), where A/, AJI and AIII are as set out in Eqn.(14), Eqn.(15) and 
Eqn.(16) respectively. 

Completing the square on Eqn.(13) leads to Eqn,(17) and substituting 
Eqn.(17) into Eqn.(5) gives Eqn.(18). Next a value for t is assigned according to 
Eqn.(19) and for r, and when u =C and u = 1 according to Eqn.(20) and 
Eqn.(21) respectively. The value of du is as shown by Eqn.(22). 

Substitution of Eqn.s(19-22) in to Eqn.(18) leads to Eqn.(23), However, 
Eqn,(24) is known and a rational approximation thereof gives the function erfix) as 
set out in Eqn,(25) with the values of r, aj, aj, aj, z^, 3^ and p as shown. 

From Eqn.(24) is finally derived the approximated tunnelling integral 
according to the method of this embodiment of the present invention and as shown 
m Eqn.(26), with the values for A/, A//, AIII, /„ ru.y(v),/(v),/"(v), A, B, C, and 
D as shown. The values of v noted above for Casel, 2 and 3 are also listed for 
Eqn.(26) as is u^^^ for v = (C+ 1)/2. Of course, u,^ is solved first in order to 
determine which of Cases 1, 2 and 3 apply. 

An additional term should be respectively added or subtracted in Eqn.(26), 
as shown in Eqn.(27), if (f„>0 & r,<0) or (r„<0 & r,>0). 

The present invention enables a commercial and practicable implementation 
of an automated simulation method for determining enhanced generation 
recombination rate due to trap-to-band tunnelling in a semiconductor device using 
the Dirac coulombic tunnelling integral. Moreover, the present invention can 
provide high levels of accuracy in the automated modelling of the enhanced 
generation-recombination rate due to trap-to-band tunnelling in a semiconductor 
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device. The method can be implemented in a two-dimensional finite element 
package. 

The present invention uses the characteristics of low, moderate and high 
field regions with a Taylor's series expansion and a reduction to one or more error 
functions. Rational approximations are applied to the error functions so as to 
provide reduction to simple exponential functions. The metod enables cancelling of 
higher order terms in the exponential functions which otherwise cause premature 
overflow errors. This enables a much wider range to be calculated, as shown in 
figure 5. 

The method of the invention enables the easy application of proper 
integration limits so as to ensure a smoth transition between low, moderate and 
high fields. The effect of removing discontinuities between the low, moderate and 
high field regions by avoiding conventional simplifications of the integration limits 
is illustrated in figure 6. 

Figure 7 illustrates actual residts of an implemention of the present 
invention for an n-channel polysilicon TFT as compared with the standard SRH 
model (no field enhancement) and the 1992 model by Hurkx et aL It is also to be 
noted that the implementation of the present invention more accurately simultes the 
leakage current at low field values (see figure 8 at Vqs = O.IV). 



wo 00/77533 



9 



PCT/GBOO/02321 



Claims 

1 . An automated simulation method for detemiining enhanced generation 
recombination rate due to trap-to-band tunnelling in a semiconductor device using 
the Dirac coulombic tunnelling integral, comprising the steps of: 

assigning the variable C to the ratio of the Poole-Frenkel barrier lowering 
energy (AE^ divided by the energy range for which tunnelling can occur (AEJ; 

assigning the value (C-h l)/2 to a variable v and performing a second order 
Taylor's series expansion of the Dirac coulombic tunnelling integral around v to 
determine a maximum value (u^) for the variable u of the integral; 

determining if the value for u^^^ is less than C, is between C and 1 or is 
more than 1; 

assigning the value of C to the variable v if is less than C; 

assigning the value of to the variable v if u^^ is between C and 1 ; 

assigning the value of 1 to the variable v if is more than 1 ; 

reducing the Taylor's series expansion of the Dirac coulombic tunnelling 
integral to an error function; 

reducing the error function to simple exponential functions by applying 
rational approximations to the error function; and 

calculating the enhanced generation recombination rate due to trap-to-band 
tunnelling in a semiconductor device using the said simple exponential functions. 

2. An automated simulation method which determines enhanced generation 
recombination rate due to trap-to-band mnnelling in a semiconductor device using 
the approximated tunnelling equation set out as equation 26 herein. 

3. An automated simulation method which determines enhanced generation 
recombination rate due to trap-to-band tunnelling in a semiconductor device using 
the approximated tunnelling equation set out as equation 27 herein. 
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4. A simulator as claimed in any one of claims 1 to 3 which determines the 
leakage current in a polysilicon Thin Film Transistor. 



5. A simulator for determining enhanced generation recombination rate due to 
trap-to-band tunnelling in a semiconductor device using the Dirac coulombic 
tunnelling integral, comprising: 

means storing a variable C having a value equal to the ratio of the Poole- 
Frenkel barrier lowering energy (AE,^ divided by the energy range for which 
tunnelling can occur (AEJ; 

means which assign the value (C-M)/2 to a variable v and perform a second 
order Taylor's series expansion of the Dirac coulombic tunnelling integral around v 
to determine a maximum value (u^^ for the variable u of the integral; 

means which determine if the value for u^ is less than C, is between C and 
1 or is more than 1; 

means which assign the value of C to the variable v if u^ is less than C; 

means which assign the value of to the variable v if u^^ is between C 

and 1; 

means which assign the value of 1 to the variable v if u^ is more than 1; 

means storing simple exponential functions derived from applying rational 
approximations to an error function obtained by reducing the Taylor's series 
expansion of the Dirac coulombic tunnelling integral; and 

. means which calculate the enhanced generation recombination rate due to 
trap-to-band tunnelling in a semiconductor device using the said simple exponential 
functions. 

6. A simulator which determines enhanced generation recombination rate due 
to trap-to-band tunnelling in a semiconductor device comprising means which 
calculate the approximated tunnelling equation set out as equation 26 herein. 

7. A simulator which determines enhanced generation recombination rate due 
to trap-to-band tunnelling in a semiconductor device comprising means which 
calculate the approximated tunnelling equation set out as equation 27 herein. 
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8. A simulator as claimed in any one of claims 5 to 7 which determines the 
leakage current in a polysilicon Thin Film Transistor. 
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